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fNJ , Abstract 

Bounds are obtained on the volume fraction in a two-dimensional body containing two 



< 



elastically isotropic materials with known bulk and shear moduli. These bounds use infor- 
mation about the average stress and strain fields, energy, determinant of the stress, and 
Q\ ' determinant of the displacement gradient, which can be determined from measurements of 

the traction and displacement at the boundary. The bounds are sharp if in each phase 
certain displacement field components are constant. The inequalities we obtain also di- 
p ! . rectly give bounds on the possible (average stress, average strain) pairs in a two-phase, 

| two-dimensional, periodic or statistically homogeneous composite. 

1 Introduction 

A fundamental problem is to determine the volume fraction occupied by an inclusion in a body, 
' or more generally the volume fraction occupied one material in a body containing two materials. 

This can usually be done by weighing the body but this may not always be practical or the 
(■ — \ densities of the two materials may be close. Then one might seek to bound the volume fraction 

' from measurements of tractions and displacements (for elasticity) or current fluxes and voltages 

(for conductivity) at the boundary of the body. If the body contains a statistically homogeneous 
or periodic composite (with microstructure much smaller than the dimensions of the body), 
then such boundary measurements can yield information about the effective tensors of the 
composite and it has long been recognized (see, for example, [281 OS E3 EH]) that bounds on 
effective tensors (which involve the volume fraction and material moduli) can be inverted to yield 
bounds on volume fractions. As shown in [21], even if the body does not contain a statistically 
homogeneous or periodic composite, but provided that the applied tractions (or boundary 
displacements), or current fluxes (or boundary voltages) are such that the fields in the body 
would be uniform were it filled with a homogeneous material, then boundary measurements can 
yield information about the effective tensor of a composite containing rescaled copies of the body 
packed to fill all space in a periodic structure. Bounds on this effective tensor yield universal 
bounds on the response of the body when such special boundary conditions are applied, which 
generalize those first obtained by Nemat-Nasser and Hori |25UH|- They can then be inverted to 
yield bounds of the volume fraction [21] , and when the volume fraction is asymptotically small 
the resulting bounds include those obtained by Capdeboscq and Vogelius [6l[7] (for conductivity) 
and Capdeboscq and Kang [8] (for elasticity) using polarizability tensor bounds [T5j[6]. Other 
bounds on the volume fraction, involving constants which are not easy to determine, were 
obtained by Kang, Seo and Sheen [13], Ikehata [12], Alessandrini and Rosset [2], Alessandrini, 
Rosset and Seo [3], and Alessandrini, Morassi and Rosset pp. 
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Given this close connection between bounding the effective tensor of a composite material 
and bounding the response of the body when these special boundary conditions are used, one 
might wonder if methods that are used to obtain bounds on effective tensors of composites could 
also be used to bound the response of the body with any boundary conditions on the fields, and 
then be inverted to bound the volume fraction. For conductivity such an approach has been 
successfully taken by Kang, Kim, and one of the authors [13], using the translation method 
of Murat and Tartar [321 [2H [33] and Lurie and Cherkaev [HI [T7] which is one of the most 
successful methods for bounding effective tensors of composites: see the books [9j HI [201 EE]- 
For a certain class of inclusion shapes (such that the field inside the inclusion is uniform for 
appropriate boundary conditions) one of the resulting volume fraction bounds gives the exact 
volume fraction. For special boundary conditions the bounds reduce to those obtained in |21j . 
and for asymptotically small volume fractions the bounds reduce to those of Capdeboscq and 
Vogelius 0E|. 

The initial goal of this paper was to use the translation method to extend [13] to elasticity, 
using measurements obtained under a single, but arbitrary, loading, and that is essentially done 
in Section [3l The attainability conditions for the resulting bounds obtained in Section [J] then 
lead us to a new method for obtaining bounds which is not based on variational principles. 
This method, the method of splitting, is described in Section [5) It correlates volume averages 
of various quantities over each phase and then projects the information onto the quantities of 
interest to obtain the desired bounds. This approach is likely to have wider applications, and in 
particular could lead to new bounds on the response of (possibly non-linear) multiphase bodies 
for many properties, not just for elasticity. 

The bounds we derive also directly give bounds on the possible (average stress, average 
strain) pairs in a two-phase, two-dimensional, composite. These bounds are the natural gen- 
eralization to elasticity of the conductivity bounds on possible (average electric field, average 
current field) pairs obtained by Raitum [291 [30] an d Tartar [M] (see also chapter 22.4 in [20 J, 
which can also be generalized to non-linear materials [221 \31\ l26j [271 [5] . Only when one of the 
phases is void, has the complete characterization of possible (average stress, average strain) 
pairs been obtained [23J. In principle bounds on the possible (average stress, average strain) 
pairs could be obtained from knowledge of the G-closure of all possible effective elasticity ten- 
sors associated with composites of the two phases mixed in prescribed proportions. However, 
this G-closure is only partly known (for a survey of results see [9l 14^ [20].) 

2 Preliminaries 

Let f2 be a smooth bounded domain in R 2 , occupied by a two dimensional elastic body. Assume 
that the body is made from two different isotropic elastic materials, characterized by their bulk 
moduli K\ and ki and shear moduli [i\ and fj,2- Denote, respectively, by a and e the stress and 
strain fields acting on the body under consideration. These fields are governed by the following 
equations. 



where u is the displacement field, S is the compliance tensor, / is the second order identity 
tensor, k is the bulk modulus (taking values k\ or K2) and fj, is the shear modulus (taking 
values Hi or 112) ■ Although two dimensional bodies do not occur in practice, this formulation 
is applicable to problems of plane stress or plane strain. 

It is assumed that one can measure the traction a-n and the displacement u at the boundary 
of Q, where n is the unit outward normal vector. From these measurements one can determine 



e = Sa = (l/2/i)<r + (1/4k - 1/4/i) (Tr <r)I, V • a = 0, 




(2.1) 
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the volume averages of certain quantities. These quantities are null-Lagrangians, functionals of 
u and/or a which can be integrated by parts, and expressed only in terms of boundary values: 
thus their Lagrangian vanishes. Introducing angular brackets to denote a volume average, i.e. 

and choosing notations so that Vu has Viti and Vu 2 as its first and second columns, rather 
than rows, the five null-Lagrangians we will work with are the average fields, 

Jan |"| Jan 

the energy, 

= lk / (*-n)-u, (2.4) 
Jan 

and the two additional null-Lagrangians 

a = (deter), 6=(detVu). (2.5) 

To express the last two quantities in terms of boundary values, it is helpful to let j% and j 2 
denote the divergence free vector fields which are the first and second columns of a, and to 
introduce the matrix 

*-(-. <26) 

for a clockwise 90° rotation. Then, as follows directly from the analysis in |13| . 



a = (ji • R±j 2 ) = T^T / qi(x)( f q 2 ), (2.7) 
Jan Jxo 

where q\ and q 2 are the fluxes q\ = ji ■ n and q 2 = j 2 ■ n, which are components of the traction 
a ■ n, xo G dtl and the last integral in (|2.7p is along the boundary 90 in the counterclockwise 
direction. Also, as follows directly from the analysis in |13j . 

b = (Vui • R±Vu 2 ) = f Ul n-R ± X7u 2 = ^- [ m^, (2.8) 

Jan |"| Jan ot 

where d/dt denotes the tangential derivative along d£l in the counterclockwise direction. 
These are not the only null-Lagrangians. For i = 1, 2 

(a ■ Vm) = — [ (a- n) Ui (2.9) 
|"| Jan 

is also a null-Lagrangian, but to simplify the analysis we refrain from considering these null- 
Lagrangians except in that linear combination which gives the energy. 

Our objective is to find inequalities (bounds) which link the values of the five null-Lagrangians 
with the volume fraction in the body and the moduli of the materials. When the body is the 
unit cell of a periodic composite material and periodic boundary conditions on the fields are 
imposed then the values of the average stress (a) and average displacement gradient (Vu) 
determine the values of the energy and a and b: 

(a ■ e) = (a) ■ (e), a = det(a), 6 = det(Vn), (2.10) 

as can be shown using Fourier analysis (see, for example, section 13.3 in [20]). Also without loss 
of generality, by making an infinitesimal global rotation if necessary, we can assume that (Vu) 
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is symmetric, in which case it can be identified with the average strain (e). Thus for composites 
the inequalities we obtain give bounds on the possible (average stress, average strain) pairs, 
and only incorporate the volume fractions of the two isotropic phases and their moduli. 

In Cartesian coordinates, a and e can be represented as 2 x 2 matrices. It is convenient to 
use the basis 





' 1 " 


1 


' 1 


1 


' 1 " 




{V2 


1 


'71 


-1 


'71 


1 


1 



so that an arbitrary symmetric matrix A 



an a\2 
0-12 0,22 



(2.11) 



is represented by v = "^( a n + 



022,011 — 022,2012). Thus, from now on we understand a and e as 3 dimensional vectors. If a 
3— dimensional vector v represents a 2 x 2 symmetric matrix A in this basis then its determinant 
is given by 

(2.12) 



det A = -v ■ Tv, 
2 



where 



1 



















(2.13) 



The null Lagrangian det a plays a pivotal role in developing the bounds of the next section. 

Recall that k and \i denote the bulk and shear, respectively, moduli of the elastic body 
under consideration. Then, in the basis £>, the tensor S can be expressed by the matrix 



r i 



S = 



- 








1 



(2.14) 



because the elastic body is isotropic. 



3 Bounds obtained by the translation method using the null- 
Lagragian det a 

Assume, in this section, that (a ■ e), do = (a), Eq = (e) and a = (deter) can be evaluated or 
estimated. For i = 1,2, let m and m be the bulk and shear, respectively, moduli of the i— th 
phase. Then, 

« = Xi«i +X2^2 and \i = xiVi (3.1) 

where 

f 1 in phase i, 
Xl \ otherwise. { ' 

Define 

fj,* = max{/xi, ^2}, k* = max{fci, K2}, (3-3) 
and fix a 6 (— ^r) to ensure that the translated tensor 

L = S-aT (3.4) 
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is positive definite. The classical complementary energy minimization principle implies 

(a ■ e) = min (a- Sa), (3-5) 

(<1> = c 
(Sa) = e 
V-a = 
(deta) = a 
a- n = a ■ n on d£l 

where the additional constraints that (a) = ao, (Sa) = £q, and (deter) = a have been added 
since we know this information about the minimizing fields. In the third constraint in f)3.5j) the 
divergence of a is understood as the divergence of the matrix that a represents. Using f|2. 12|) 
we then have 

(a ■ e) — 2aa = min (a • La). (3-6) 

(c) = a 
(Sa) = e 
V -a = 
(det a) = a 
a - n = a ■ n on d£l 

Dropping the last three constraints in the minimum above and defining 

eo = e - aTa , (3.7) 

gives 

(a • e) — 2aa > min (a -La). (3-8) 
(c) = a 
(La) = e 

This minimum can be found using the Lagrange multiplier method. In fact, if a denotes 
the minimizer of the right-hand side of (|3.8p . then there exist two constant vectors Ai and A2 
such that 

2La = X 1 + LX 2 . (3.9) 

Denoting by (a)i and (17)2 the average of a on phase 1 and phase 2, respectively, from (|3.9p . 
we have 

2Lx(a)i = A1 + L1A2, 

2L 2 (a) 2 = A1 + L2A2, (3.10) 

which gives 

Ai = 2(Li - L 2 )- 1 L 1 L 2 ((a) 2 - (a)i), 

A 2 = 2(L 1 -L 2 y l (L 1 (a) 1 -L 2 (a) 2 ). (3.11) 
On the other hand, we can see that 

(cr)i = ~r(Li — L 2 )~ 1 (eo — L 2 ao), 
Ii 

(ff>2 = -^(Li-L 2 )- 1 (eo-Lia ), (3.12) 

J2 



by solving the system 



co = /i(c)i + f 2 (a) 2 

e = fiL 1 (a) 1 + f 2 L 2 (a) 2 . (3.13) 
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Thus, we have 

Ai = -^-{L l -L 2 )- 2 L l L 2 {e G -{L)a Q ), 

A 2 = 7 ^( J Ci-L 2 )- 2 [(/ 2 L 1 + / 1 L 2 ) eo -L 1 L 2 a ]. (3.14) 

Jl/2 

For simplicity in calculations, assume first that a G (— , 2^-) so that L is invertible. Later 
we will consider the limiting case where L is singular, but still positive semi-definite. It follows 
from {331 tnat 

2e = Ai + (L)A 2 , 

2^ = (L- 1 )Ai + A 2 . (3.15) 

Moreover, we have 

{a -La) = i((L- 1 Ai + A 2 )-(Ai + LA 2 )) 

= i[(L- 1 )A 1 -Ai + 2Ai-A 2 + (L)A 2 -A 2 ]. (3.16) 

On the other hand, it follows from (|3.15p that 

Ai • A 2 = 2s • A 2 - (L)A 2 • A 2 = 2a • Ai - (L -1 )Ai • Ai. (3.17) 
Thus, it follows that 

(a ■ La) = - [a ■ Ai + A 2 • e ] 

= i[(2eo-(L)A 2 )-ao + A 2 -e ] 
(L — L )~ 2 

= e ■ * + 1 1 , , V [eo - (L)a ] • [(/ 2 Li + / x L 2 )e - £iL 2 <tq]- (3.18) 
/1/2 



This, together with (|3.7() and (j3.8|) . gives us the bound 

(L — L )~ 2 

(a-e)-a -e - 2aa + 2a det cr > — 1 , , 2 [e - (i)cro] ■ [(/ 2 ^i + fiL 2 )e - LiL 2 a ] . (3. 19) 

Jl/2 

By taking limits, we can see that (|3.19p is valid for a not only in ( — 2^*' 2 «*) but a ^ so a ^ 
— 2^* and 2^*. Although (|3. 19|) looks like a quadratic inequality with respect to a (because 
of the definitions of eo and L in (|3.7p and (|3.4p respectively), we can show that it is linear 
by expanding (fiL 2 + f 2 L\)eo — L\L 2 o~q and seeing that the coefficient of a 2 is 0. Hence, the 
bound in ()3.19|) improves or gets worse (depending on the data given) as a tends to — or 
2^ with the optimum value for a occurring at one of the two limits. The arguments above can 
be summarized as 

Theorem 3.1 The following bound 
{a ■ e) - (j • e 

> min {t^(5i - S 2 )- 2 (e - (S)a ) ■ {{hL 2 , + f 2 L u )e - L u L 2 ,a ) - ^f^, (3.20) 
- <S 2 )- 2 (e - (S)a ) ■ {(f x L 2 * + / 2 Li*)e - L^L 2 *a ) + , 

where L; = LA n _ 1 and Li* = LA n _ 1 , i = 1,2, holds. 
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4 The attainability condition for (13.201) 



In this section, we find conditions for the field a to be such that the equality for (|3.20p is 
attained. In other words, fixing a G [— ^ , ^*], we find conditions on a such that 

(a ■ La) = (a -La), (4.1) 

where a is a minimizer of the right-hand side in (|3.8p which can be found from (|3.9p and (|3. 1 1 j) . 
We have the following lemma. 

Lemma 4.1 J^. 1\ ) holds if, and only if, 

La = La = -7-7- (I-i - L 2 )~ 1 [f 2 xiL 1 (e - L 2 a ) - fiX2L 2 (e - Li<x )]. (4.2) 
Jih 

Proof: Since L is self-adjoint, the "only if direction is not hard to see. We thus only prove 
the "if direction. Define the functional J by 

J(a) = {a-La), (4.3) 

for all vector valued functions a satisfying 

(ff) = o-o, {La) = e . (4.4) 

Since L is semi-positive definite, J is convex. It is not hard to see that for all t £ [0,1], 
a = ta + (1 — t)a satisfies (|4.4p . Using (|4.ip and the fact that a and a are both minimizers of 
J, we have 

J{a) < J(ta + (1 - t)a) < tJ{a) + (1 - t)J{a) = J(a) (4.5) 
for all t G [0, 1] . It follows that 

j t J(ta + (l-t)a) = 0, (4.6) 

or equivalently 

(L(ta + (1 - t)&) ■ (a - a)) = 0. (4.7) 
Letting t in (|4.7p be 1 and 0, respectively, gives 

(La ■ (a - a)) = (4.8) 

and 

(La ■ (a - a)) = (4.9) 
From the difference of these two equations, 

((a - a) ■ L(a - a)) = 0, (4.10) 

we obtain the first equation of (|4.2p . 

In order to see the second equation of (14. 2j) . we calculate La on each phase. On phase 1, 

La = Ai + L\X 2 

= 2e - (L)X 2 + L x \ 2 

= 2e + f 2 (L 1 -L 2 )X 2 
2 

= 2e + -r{L\ - ^ 2 ) _1 [(/2^i + fiL 2 )e - L x L 2 a Q \ 
Ji 

2 

= —(L 1 -L 2 ) [/i(Li - L 2 )e + (/2L1 + fiL 2 )e - LiL 2 a Q \ 
Ii 

= ^{L x -L 2 )- l L x {e^-L 2 a Q ). (4.11) 
/1 
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Similarly, on phase 2, 

La = ~{L X - L 2 )- l L 2 (e Q - L l a Q ). (4.12) 

J2 

The lemma follows. □ 
Lemma 14.11 implies that if bound (|3.20[) is attained then there are two constant vectors D 
and E such that 

La = La = xiD + X 2E. (4.13) 

We next show that (I4.13P is the attainability condition for (I3.20P , that we are looking for in this 
section. In fact, denoting by (a)i and (a) 2 the average of a on phase 1 and phase 2, respectively, 
as in the previous section, we have 

D = L 1 {a) 1 , E = L 2 (a) 2 . (4.14) 

Plugging (I3.12j) into (|4.14p and then the resulting values of D and E into (I4.13p . we obtain 

(32}. 

We have proved the theorem. 
Theorem 4.2 The bound in \3.20) becomes equality if, and only if, 

L,a = xiD + X 2E (4.15) 

where D and E are two constant vectors and L* is either LI i or L\ n _ i . Moreover, 

1. if fi* = [i\ 7^ \x 2 then 
{a ■ e) - a ■ e 

is equivalent to 



L\ a= —L.o = XlA. + X2E (4.17) 

with D* = (rfi,0, 0) and this holds if and only if the field a is constant on phase 2, and 
has constant first component (bulk part) on phase 1; 

2. if k* = K\ k 2 then 

(a ■ e) - a ■ e 

jk- 2 (Si - S 2 )- 2 (e - (S)a ) ■ {(hL 2 , + f 2 L x *)e - L v ,L 2 ,a a ) + 



~ fif: 

is equivalent to 



L\ a= ^o = X iD* + X2E (4.19) 



with D t . = (0, d 2 , d^) and this holds if and only if the field a is constant on phase 2, and 
has constant second and third components (shear part) on phase 1;. 
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5 The method of splitting 



We introduce here another approach, not based on variational principles, using which we can 
deduce the previous bounds. The main idea of the method is to split the domain into its 
phases, to correlate averages over each phase, and then to project the information to obtain the 
desired bound. Moreover, this method allows us to add one more datum: the null-Lagrangian 
b = (det Vii). In other words, the known quantities are 



E = (a-Vu), a = {a), (Vu), a = (deter), 6=(detVu) 
Consider the four dimensional space of 2 x 2 matrices with the basis 





1 " 


1 


' 1 " 


1 


' 1 





1 


' 1 " 






-1 


'71 


1 


'71 





-1 


'71 


1 


! 



so that 
and 



a = (0,0-1,0-2,0-3), Vn = (F ,£i,£2,£3) 
\(o{ - a\ - al), c^b- i(F ) 2 > b - i(F 2 ) = 1 - 



2 N u ' 2 

Note that the inequality in (|5.4p is attained if, and only if, F$ is constant everywhere. 
Defining (b for bulk and s for shear) 



E2s 



we have 



= (xio-iei) = 2ki(xi£i), 

= (X2^i£i) = 2/c 2 (X2£i), 

= (xi(°2£2 + cr 3 e 3 )) = 2ui(xi(ei + £§))» 

= (X2(<72£2 + a 3 e 3 )) = 2u 2 (X2(ei + 

E = Eu + -Ei s + E2b + E2 S , 
a = K\Eib + K2E 2 b — Ml -^ls — M2-E 2s , 

E\b E-ib Eis E2s 

c > ! . 

4k 2 4k 2 4ui 4/i 2 



On the other hand, it follows from the identities 



that 



1 



2(ac 2 — Kl 

1 

2(ji 2 -m) 



-(2k 2 (£i) - (0-1)), (X2£l> 
(2^2 (ej) ~ Wj)), (Xisj) 



1 



2{Kl - K2) 
1 

2(u 2 - ui) 



(5.1) 

(5.2) 

(5.3) 
(5.4) 



(5.5) 

(5.6) 
(5.7) 

(5.8) 



(Xl £ i) + (X2£i), 2 = 1,2,3, 

2/ci(xi£i> + 2/t 2 (x 2 e 2 ), (<Tj) = 2ui(xi£j> + 2u 2 (x2£j>, j = 2, 3, (5.9) 



(2/ci(ei) - (o-i)), 

(2u 2 (e,) - (o,-)), j = ^5310) 



So, these quantities are known. 
Note that for i = 1, 2, 3, 



(Xief) " T (xi£i} 2 
h 



((Xi£i 



^(Xi^)) 2 ) > 0, 
Ji 



(5.11) 
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with equality when is constant on phase 1. Similarly, 



(X2£?> - ^(X2^) 2 > 0, i = 1,2,3, (5.12) 

/2 



Ai6 = 2ki(xi£i) 2 , 

A 26 = 2fi; 2 (X2ei) 2 , 

A u = 2 f i 1 ({ X ie 2 ) 2 + {xie 3 ) 2 ), 

A 2s = 2^ 2 ((X2£2) 2 + (X2£ 3 ) 2 ), 



with equality when e$ is constant on phase 2. Therefore, defining the known quantities 

(5.13) 

■Ms — ^PU\A1 C 2/ T- \Al c a/ 7) 

4 2s = 2/x 2 ((x 2 £ 2 ) 2 + (x2£3> 2 ), 
we have from (|5. 11 j) and (|5.12p 

A . A.. 

(5.14) 
(5.15) 



Eu > 


A lb 
h ' 


E 2b > 


A 2b 
h ' 


Eis > 


h ' 


E 2s > 


A 2s 
h 



Also, E\ b and E 2b can be eliminated using (|5.6p and (|5.7p 



^16 = ^7^( a ~ + («2 + Hl) E ls + (« 2 + fJ-2)E 2s ), 

E 2 b = (a - Ki-E + («! + m)E ls + («! + n 2 )E 2s ). 



So, (|5.14p gets replaced by 



and the inequality in (|5.8j) gets replaced by 

C - 4ki(<ci-k 2 ) ( Q ~ + ( K2 + Ml)- 5 ^ + ( K 2 + ^2)E 2s ) 

+T^k^D (° - KlE + ( Kl + + ( Kl + - if - if ■ 

Now we have 

1 1 1 

+ 



(5.16) 



1 Aih 

(a - k 2 E + (k 2 + + (k 2 + /U 2 ) J B 2s ) > -rS (5.17) 

Kl - K 2 /l 

and ^ 

(a - k x E + (ki + ^1)^1, + («i + M2)^ 2s ) > (5.18) 

k 2 - ki h 



(5.19) 



ki(ki-k 2 ) k 2 (k 2 -ki) KlK 2 

—K 2 _ Kl Kl + K 2 

- K 2 ) K 2 (K 2 — K\) K\K 2 

K 2 +Hl Kl+/il 1 _ HiK 2 (K2 + Hi) — HiK\{ki + Hi) — kik 2 (ki — K 2 ) 



ki(ki — k 2 ) k 2 (k 2 — ki) hi kik 2 hi{ki — K2) 

-jmXki + K 2 ) -Hl~ K l K 2 
Kl«2A*l 

(ki + hi)( k 2 + A^i) 
K1K2M1 

K1+M2 | K 1 + H2 1 _ («1 + M2)(^2 + Mg) /g 



Ki(ki — K 2 ) K 2 (k 2 — «l) /Z2 KiK 2 H 2 
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So, (|5.19p becomes 



Ex s (K\ + ^l)(«2 + Ml) . ^2s(«l + M2)(«2 + M2) 



/il 



+ 



/i 2 



> — AK1K2C. 



(5.21) 



So, (|5,15|) . (|5.17|) . (|5.18|) and (|5.21|) give us 5 inequalities on the pair (E\ s , E2 S ). To project out 
the information about the unknowns (Ei s , E 2s ) we observe that the volume fraction /1 = 1 — /2 
must be such that these inequalities define a feasible region in the (E\ s ,E 2s ) plane (see e.g. 
Figure [1]). Note that when either f\ or f 2 goes to 0, (|5. 15[) cannot be satisfied. In this case, 

E 2s 

The line defined by (15.21h 




The line defined by (|5T7j) 
The line defined by (|5.18p 



Eis 

Figure 1: The region of possible (Ei s ,E2 S ), assuming that k% > k 2 . 

the feasible region is empty. Thus, in the generic case, at a bound on f±, i.e. at the limiting 
value of fi, the feasible region shrinks to a point. In other words, in the generic case, 3 of 
the inequalities will be satisfied as equalities and the remaining 2 as inequalities: the picture 
looks like, e.g. Figure [2j Note that the 3 inequalities which are satisfied as equalities must have 

E 2s 



The line defined by (|5.21|) 




he line defined by (|5.17p 
The line defined by (Qgj) 



Eu 

Figure 2: The feasible region shrinking to a point, assuming that K\ > K 2 . 
"outward normals" ni,n 2 ,ns such that 

a\n\ + a 2 ^2 + 03^3 = 



(5.22) 



for some a±,a 2 ,a3 > 0. The previous approach in Section [3] only took the four inequalities in 
(I5.15p . (I5.17P and (15.180 into account. The limiting value of fi corresponds to the case where 3 
of these were satisfied as equalities and the fourth as an inequality. Thus, one sees immediately 
that the attainability condition is where the field is constant in one phase and either the bulk 
part or shear part is constant in the other phase. With the additional inequality (|5.2ip . there 
will be 10 cases to consider (number of ways to pick 3 equations from 5), 4 considered above 
and 6 new. However, not all of these 10 cases satisfy (|5.22p ; e.g. if K\ > k 2 , the triplet of 
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equations that correspond to (|5. 15|) and (|5.2ip do not satisfy (|5.22p since if E\ s and E<i s are 
both large and positive, all three equations are satisfied; i.e. the feasible region given by the 
three equations cannot be a point, but is instead an open region. The question arises as to 
which of these 6 new cases satisfy (|5.22p . 

It is obvious that the outward normals of the regions that satisfy the two inequalities in 
(|5.15p are, respectively, 

i^i = (—1,0) (5.23) 

and 

V2 = (0,-1). (5.24) 

Without loss of generality, by relabeling the phases if necessary, assume that k\ > K2- In this 
case, an outward normal of the region that satisfies (|5.17p is 

^3 = (~«2 - Ml)-«2 - A*2), (5.25) 

an outward normal of the region that satisfies (|5.18p is 

i>4 = («i + + M2), (5.26) 

and, finally, an outward normal of the region that satisfies (|5.2ip is 

v = ( («1+Ml)(«2+Ml) (^i +^ 2 )(K 2 +/i 2 ) \ , 5 27 ^ 

We now consider each triplet of normals that include v§ to see when they satisfy (|5.22|) . As 
said, (^1,^2,^5) does not satisfy (|5.22p . If (a,f3) solves 

u 5 = av x + pV 3 , (5.28) 

then 

P = > 0. (5.29) 

M2 

Thus, the triplet [y\, U3, u§) never satisfies (|5.22p . Similarly, neither does the triplet (f 2 , ^3, 1/5). 
Next, solving the equation 

v 5 = av\ + pV 4 , (5.30) 

gives 

a = _«aO*i-M2)(«i + Mi) p = _"* + l» m (5.31) 
M1M2 M2 
So, the triplet (1/1, 1/4, v§) satisfies (|5.22p if, and only if, \i\ > /x 2 . By a similar technique, solving 

1/5 = cw 2 + pV 4 , (5.32) 

we have 

a = _"s(«i+^)(M2-Mi) /3 = - K2 + ^. (5.33) 
M1M2 Mi 
Therefore, the triplet (^,^4,^5) satisfies (|5.22|) if, and only if, \x\ < n%. Finally, we consider 
the triplet (^3,1^4,^5). The unique solution of 

1/5 = av^ + /3z/4 (5.34) 

is 

«2(«i +Mi)(«i + M2) ^ n o «l(«2 + MlK^2 + M2) ^ n 
a = 7 r < U, p = 7 r < U. (5.35) 

(«l-«2)MlM2 (Kl-«2jMlM2 

Hence, this triplet always satisfies (|5.22p . 

We next find the desired bounds. We can do so by finding an appropriate linear combination 
of the equations in each triplet involving (|5,2ip that satisfies (|5,22p to obtain 
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Theorem 5.1 Assuming K\ > k 2) and recalling the definitions A5.13\) and \5.l 0\) of An,, A2b, 
A\ s , and A2 S one has the bound 

KiK 2 {a + E(ni + IX 2 )) , A Xh K 2 (Kl + Hi)(Kl + Hi) , ^2&«l(«2 + A*l)(«2 + H2) 

4CK1K 2 > 1 J 1 J , 

fJ-1^2 71^1^2 hV\V>2 

(5.36) 

with equality when Fq ( the component of the antisymmetric part of Vu ) is constant everywhere, 
and the bulk component e\ is constant in phase 1 and constant in phase2. If additionally m > \i<i 
then one has 

K 2 {a + E(-ki + fi 2 )) , A 1s k 2 (k 1 - H2) A 2b (n 1 - k 2 ){k 2 + /u 2 ) 

ACK1K2 > 1 7 1 7 ; (5.37) 

flfilfJ-2 J2V2 

with equality when Fq is constant everywhere, the shear components £2 and £3 are constant 
in phase 1, and the bulk component e% is constant in phase 2. Alternatively if Hi < ^2 (and 
ki > K2) then one has 

K 2 {a + E(-K 1 + Mi)) A 2s K 2 (pk2 - An)(«i + H2) , A 2b (K 1 - k 2 ){k 2 + m) 00 . 

4ckik 2 > 1 7 1 7 > (5.38) 

A*l 11^2 J2H1 

with equality when Fq is constant everywhere, and both bulk and shear components £\, e 2 and 
£3 are constant in phase 2. 

Proof: As suggested by (EMI) , we can add (f5T2TD to d5T7|> . multiplied by K ^ Kl + + 

and (|5.18p . multiplied by i( K 2 + Mi)( 2 + ^2) ^ eliminate both E\ s and £" 2s . Doing so, we 

obtain (15.36p . Similarly in the case that |Ux > /^2j we can use the first inequality in (15,15p . fj5 .181) 
and (|5.2ip to deduce (|5,37p . Finally, if fi\ < ^ 2 then the second inequality in (|5.15p . (|5.18p 
and (15.2ip yield (I5.38p . The attainability conditions follow directly from the conditions under 
which the inequalities in (|5.4p and (|5.1ip are satisfied as equalities. □ 
It is an open question as to whether the bounds of Theorem 15.11 could be obtained using the 
translation method. Although we suspect they could, the method of splitting has the advantage 
of immediately providing the attainability conditions. 

Remark 1 The final set of bounds are the intersection of the inequalities provided by Theorem 
\3.1\ and Theorem \5.1{ Each of these inequalities can be multiplied by /i/ 2 to yield quadratic 
inequalities in fi = 1 — f 2 which may be easily solved to give the maximum interval of f\ 
compatible with all the inequalities. If one is interested in bounds on (average stress, average 
strain) pairs in composites then one should make the substitutions \2.10) in these inequalities. 
It remains to be investigated whether the resulting bounds provide a complete characterization 
of the possible (average stress, average strain) pairs in composites of two phases mixed in a 
given proportion, or whether there are some missing bounds. 
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